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1. Introduction 

The hard dynamics of processes involving hadrons is nowadays remarkably well de- 
scribed by QCD predictions. An important role has been played in this achievement 
by the ability to compute the relevant reactions to next-to-leading order (NLO) ac- 
curacy, a mandatory step in view of the large value of the coupling constant a s . 
Although fairly successful in their phenomenological applications, NLO predictions 
are affected by uncertainties that can be of the order of a few tens of percent, typ- 
ically estimated by varying the renormalization and factorization scales. This is a 
consequence of the fact that, in some cases, NLO corrections are numerically as im- 
portant as the leading order (LO) contributions. By increasing the accuracy of the 
perturbative predictions, through the computation of the next-to-next-to-leading or- 
der (NNLO) contributions, one would certainly reduce the size of the uncertainties, 
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and obtain a firmer estimate of the rates [1]. Such a task, however, implies finding 
the solution of a few highly non trivial technical problems. Recently, two major 
bottlenecks have been cleared: two-loop functions with up to four legs have been 
computed (with zero or one massive leg) [2], and so have the three-loop Altarelli- 
Parisi kernels [3], which opens the way to exact NNLO PDF fits. 

The present situation is thus fairly similar to that of the early 80 's. Back then, the 
computation of the two- loop Altarelli-Parisi kernels [4], and the ability to compute 
all the tree-level and one-loop functions for specific processes, left open the problem 
of achieving explicitly the cancellation of soft and collinear singularities, as predicted 
by the Kinoshita-Lee-Nauenberg (KLN) theorem for inclusive, infrared-safe observ- 
ables. Early approaches pioneered the subtraction [5] and phase-space slicing [6] 
techniques, for computing analytically the divergent part of the real corrections, in 
the context of the prediction for a given observable; a different observable required a 
novel computation. Later, it was realized that the KLN cancellation could be proven 
in an observable- and process-independent manner, which allows one to predict any 
observable (for which the relevant matrix elements can be computed) to NLO accu- 
racy without actually using the observable definition in the intermediate steps of the 
computation [7]-[10]. Apart from being very flexible, these universal methods have 
the virtue of clarifying the fundamental structure of the soft and collinear regimes 
of QCD. 

Observable- and process-specific calculations have certain advantages over uni- 
versal formalisms; the possibility of exploiting the observable definition and the kine- 
matics of a given reaction in the intermediate steps of the computation generally 
leads to more compact expressions to integrate. This is evident if we consider the 
fact that the first NNLO results for total rates [11, 12] pre-date the universal NLO 
formalisms by some years. More recently, the first calculation of a rapidity distri- 
bution at NNLO has been performed using these techniques, for single vector boson 
hadroproduction [13, 14]. 

In the last two years a fair amount of work has been carried out in the context 
of process-specific calculations. The structure of NNLO infrared-singular contri- 
butions necessary to implement the subtraction method has been discussed in the 
case of e + e" — > 2 jets [15], and the C% contribution to e + e~ — > 3 jets has been 
computed [16]. Moreover, an observable-independent method [17], based on sector 
decomposition [18], has been recently proposed. This approach allows one to handle 
and cancel infrared singular contributions appearing in the intermediate steps of an 
NNLO calculation in a fully automatic (numerical) way, and thus significantly differs 
from the semi-analytical approaches of refs. [7]-[10]. This method has been applied 
to the computations of e + e~ — > 2 jets [19] and Higgs hadroproduction [20] cross 
sections. 

One may argue that, since the number of two-loop amplitudes computed so far is 
limited, process-specific (but observable-independent) computations such as those of 
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refs. [16, 19, 20] are all what is needed for phenomenology for several years to come. 
Although this is a legitimate claim, we believe that universal formalisms are interest- 
ing in themselves, and that their achievement should be pursued in parallel to, and 
independently of, that of process-specific computations. Work in this direction is 
currently being performed by various groups, and partial results are becoming avail- 
able [21]-[28]. We stress that the formulation of universal subtraction formalisms at 
NNLO would pave the way to the matching between fixed-order computations and 
parton-shower simulations, similarly to what recently done at NLO [29, 30], thus 
resulting in predictions with a much broader range of applicability for phenomeno- 
logical studies. 

The purpose of the present paper is to propose a general framework for the 
implementation of an observable- and process-independent subtraction method. A 
complete formalism would require the construction of all of the NNLO counterterms 
necessary to cancel soft and collinear singularities in the intermediate steps of the 
calculation, as well as their integration over the corresponding phase spaces. Most of 
the recent work on the subject concentrated on the former aspect, without performing 
the analytical integration of the kernels. In this paper we take a different attitude: we 
propose a subtraction formula, and we give a general prescription for the construction 
of all of the counterterms. However, we do not construct most of them explicitly 
here, since we limit ourselves to considering only those relevant to the C F T R part 
of the dijet cross section in e + e~ collisions. On the other hand, we specifically 
address the problem of their integration. We propose general formulae for the phase 
spaces necessary for the integration of the subtraction kernels, and we use them to 
integrate the counterterms mentioned above. In this way, we directly prove that the 
subtraction formula we propose does allow us to cancel the singularities at least in 
the simple case of the C F T R contribution to the e + e~ — > 2 jets cross section, and we 
construct a numerical code with which we recover the known results for dijet and 
total rates. Since the subtraction formula, the counterterms, and the phase-space 
measures are introduced in a way which is fully independent of the hard process 
considered, this result gives us confidence that the framework we propose here is 
general enough to lead to the formulation of a complete subtraction formalism, whose 
explicit construction is, however, beyond the scope of the present paper. 

The paper is organized as follows: in sect. 2 we review the strategy adopted by 
NLO universal subtraction formalisms, and we formulate it in language suited to its 
extension to NNLO. Such an extension is discussed in sect. 3, and our main NNLO 
subtraction formula is introduced there. This subtraction procedure is shown to work 
in a simple case in sect. 4. We present a short discussion in sect. 5, and give our 
conclusions in sect. 6. A few useful formulae are collected in appendices A and B. 
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2. Anatomy of subtraction at NLO 



Let us denote by r any real (tree-level) matrix element squared contributing to 
the NLO correction of a given process, possibly times a measurement function that 
defines an infrared-safe observable, and by d$ the phase space. For example, when 
considering two-jet production in e + e~ collisions, r is the product of the 7* — > qqg 
matrix element squared times the S functions that serve to define the jets within a 
given jet-finding algorithm, and rf$ is the three-body phase-space for the final-state 
partons q, q, and g. As is well known, the integral 



is in general impossible to compute analytically. In the context of the subtraction 
method, one rewrites eq. (2.1) as follows: 



The quantities f and d& are completely arbitrary, except for the fact that fd$ must 
fulfill two conditions: the first integral in eq. (2.2) must be finite, and the second in- 
tegral must be calculable analytically; when this happens, fd& is called a subtraction 
counterterm. The former condition guarantees that the divergences of frd<& are the 
same as those of JrdQ. Thus, by computing JfdQ one is able to cancel explicitly, 
without any numerical inaccuracies, the divergences of the one-loop contribution, as 
stated by the KLN theorem; when using dimensional regularization with 



the divergences will appear as poles l/e fc , with k — 1, 2. As far as the first integral in 
eq. (2.2) is concerned, its computation is still unfeasible analytically. However, being 
finite, one is allowed to remove the regularization, by letting e — * 0, and to compute 
it numerically. 

In the context of the subtraction method, the computation of an observable to 
NLO accuracy therefore amounts to finding suitable forms for f and that fulfill 
the conditions given above. Clearly, it is the asymptotic behaviour of r in the soft 
and collinear configurations that dictates the form of f. The major difficulty here 
is that soft and collinear singularities overlap, and f should not oversubtract them. 
To show how to construct f systematically, let us introduce a few notations. We call 
singular limits those configurations of four-momenta which may lead to divergences 
of the real matrix elements. Whether the matrix elements actually diverge in a given 
limit depends on the identities of the partons whose momenta are involved in the 
limit; for example, no singularities are generated when a quark is soft, or when a 
quark and an antiquark of different flavours are collinear. We call singular (partonic) 
configurations the singular limits associated with given sets of partons, and we denote 
them as follows: 




(2.1) 





d = 4 - 2e 



(2.3) 
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• Soft: one of the partons has vanishing energy. 

S>: P? - 0. (2.4) 



• Collinear: two partons have parallel three-momenta. 

dj : pi || pj. (2.5) 



• Soft-collinear: two partons have parallel three-momenta, and one of them is 
also soft. 

SC ir . (p?-0,a||ft-)- (2-6) 

The singular limits will be denoted by S, C, and SC, i.e. by simply removing 
the parton labels in eqs. (2.4)-(2.6). In practice, in what follows parton labels will 
be often understood, and singular limits will collectively indicate the corresponding 
partonic configurations. 

For any function /, which depends on a collection of four-momenta, we introduce 
the following operator: 

S NLO (f) = f - E f(Si) - E fidj) - E fiSQj), (2.7) 

i i<j ij 

where f(L a ), L a = S a ,C a , SC a is the asymptotic behaviour of / in the singular 
partonic configuration L a . Notice that the role of indices i and j is symmetric in the 
collinear limit, but is not symmetric in the soft-collinear one (see eqs. (2.5) and (2.6)), 
and this is reflected in the sums appearing in the third and fourth terms on the r.h.s. 
of eq. (2.7). The asymptotic behaviour f(L a ) is always defined up to non-singular 
terms; however, what follows is independent of the definition adopted for such terms. 
Although all of the singularities of / are subtracted on the r.h.s. of eq. (2.7), S NLO (f) 
is not finite, owing to the overlap of the divergences. To get rid of this overlap, we 
introduce a set of formal rules, that we call the £- prescription: 

1. Apply 

&NLO 

to /, getting f + Eaj£- 

2. Apply S NLO to every obtained in this way and substitute the result in the 
previous expression, getting / + £ Q1 /£) + J2 a2 fif- 

3. Iterate the procedure until SnloU^) = V«„, at fixed n. 

4. Define /=-EiEa i /iJ- 

In order to show explicitly how the ^-prescription works, let us apply it step by step. 
After the first iteration, we find eq. (2.7). With the second iteration, we need to 
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compute 



S NLO (f(Si)) = f(Si) - E f(Si © Si) - E f(Si © C jk ) - £ f(Si © 5C jfc ), (2.8) 

j^i j<k jk 

W/(Ctf)) = fidj) - E f(dj ®S k )- E f(dj © c«) 

k k<l 

~ E WjtBSda), (2.9) 

£ NLO {f{Sdj)) = f(Sdj) - E f(Sdj © 5 fc ) - E f(Sdj © Cm) 

fc fc</ 

- E HSdjQSCu), (2.10) 

where we denoted by f(L ai ®L a2 ) the asymptotic behaviour of the function /(L Ql ) in 
the singular partonic configuration L a2 . Although in general the operation © is non 
commutative, we shall soon encounter examples in which /(L ai ©L a2 ) = f(L a2 ®L ai ). 
Notice that we included all possible singular parton configurations in eqs. (2.8)- 
(2.10), except for the redundant ones - an example of which would be the case k — i 
in the second term on the r.h.s. of eq. (2.10). 

We now have to take into account the fact that we are performing a computation 
to NLO accuracy. Thus, the definition of an observable will eventually be encoun- 
tered (for example, embedded in a measurement function), which will kill all matrix 
element singularities associated with a partonic configuration that cannot contribute 
to the observable definition at NLO. For example, the limit in which two partons are 
soft is relevant only to beyond-NLO results, and this allows us to set f(Si © Sj) = 
V{?,j}; analogously 

J2f( C ^® S k) = f(C l3 ®S t ), (2.11) 
k 

since the case k ^ i would result into two unresolved partons, which is again a 
configuration that cannot contribute to NLO. In general, it is easy to realize that 
the operator £ NLO is equivalent to the identity when it acts on the terms generated 
in the second iteration of the ^-prescription (i.e., the terms with negative signs in 
eqs. (2.8)-(2.10)), and thus that, according to the condition in item 3 above, the 
^-prescription requires at most two iterations at NLO (it should be stressed that 
this would not be true in the case of an infrared- unsafe observable, which would lead 
to an infinite number of iterations). It follows that 

/ = E m) - E m © dj) + E fiPa) - E © + E /(sc y ). (2.12) 

i ij i<j ij ij 
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This expression can be further simplified by observing that the soft and the collinear 
limits commute. This allows one to write 

f(S t © dj) = fidj © S t ) = /(SCij), (2.13) 

and therefore 

/ = E f(Si) + E " E / W,-). (2-14) 

i i<j ij 

Let us now identify / with rrf$. Eq. (2.14) implies that the term added and sub- 
tracted in eq. (2.2) reads 

fd$ = JX^^) +Y,r(C ij )d$(C ij ) -Y^riSC l3 )d<$>(SC l3 ). (2.15) 

i i<j ij 

All that is needed for the construction of a subtraction counterterm at NLO is 
eq. (2.15), and the definition of the rules for the computations of r(L) and d$(L). In 
other words, all subtraction procedures at NLO are implementation of eq. (2.15), i.e. 
of the ^-prescription, within a computation scheme for the asymptotic behaviours of 
the matrix elements and the phase spaces. Although not strictly necessary in princi- 
ple, it is always convenient to adopt the same procedure for the computation of r(L) 
and of d<&(L); most conveniently, this is done by first choosing a parametrization for 
the phase space, and by eventually using it to obtain r(L). 




Figure 1: Configurations contributing to soft (left panel) and collinear (right panel) limits 
of squared amplitudes. 

In order to be more specific, we shall consider two explicit constructions of sub- 
traction counterterms, namely those of ref . [8] and of ref . [9] . Since the subtractions 
always need to be performed at the level of amplitudes squared, the relevant diagrams 
(in a physical gauge) are those depicted in fig. 1, for the soft (left panel) and collinear 
(right panel) limits respectively. According to the notations introduced before and 
the parton labelings that appear in the diagrams, we denote the corresponding sin- 
gular partonic configurations by Si and C^. As is well known, the leading (singular) 
behaviours of the real matrix elements squared will be given by the following factors 
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(see e.g. ref. [31]) 



Pj-Pk 



(2.16) 




Pi-PjPi-Pk 
1 



(2.17) 



Pi-Pj ' 



with the parton k playing no role in the collinear limit. In order to integrate the 
subtraction counterterms analytically, a phase space parametrization must be chosen 
such that the leading divergences displayed in eqs. (2.16) and (2.17) have as trivial as 
possible a dependence upon the integration variables. Furthermore, the eikonal and 
collinear factors of eqs. (2.16) and (2.17) have manifestly overlapping divergences; 
thus, a matching treatment of the two, while not strictly necessary, would allow an 
easy identification of the overlapping contributions. 

In ref. [8] the phase space is first decomposed in a manner which is largely arbi- 
trary, but such that in each of the resulting regions only one soft and one collinear 
singularity at most can arise (i.e., the other singularities are damped by the 9 func- 
tions which are used to achieve the partition); the two may also occur simultaneously. 
Thus, each region is identified by a pair of parton indices - say, % and j - and no 
singularity other than 1/pi-pj can occur in that region. This implies that the eikonal 
factor in eq. (2.16) will not contribute a divergence to the region above when pi || pfc. 



integration variables, but is irrelevant in the treatment of the singularities. Clearly, 
this parametrization is suggested by a collinear-like configuration, but thanks to the 
partition of the phase space it also allows a straightforward integration over soft 
singularities. Graphically, this is equivalent to squaring the parts of the diagrams in 
fig. 1 which lay to the left of the cut. The contribution due to the part to the right 
of the cut in the diagram on the left panel is clearly recovered once the sum over 
parton labels is carried out, since the role of indices j and k is fully symmetric. Once 
the exact parametrization of eq. (2.18) is fixed, ref. [8] proceeds by defining 




In turn, this makes possible to choose a parametriza- 
tion of the phase space, based on exact factorization 
formulae (see app. B), in which a pseudo-parton ij 
(thick line in fig. 2) branches into on-shell partons i 
and j (narrow lines); in other words, the phase space 
of eq. (2.1) is written as follows: 



d$ = d$ n _i(i7')d$ 2 (i,j), 



(2.18) 



Figure 2: Splitting topol- 
ogy contributing to NLO 
computations. 



where is the two-body phase space, times the mea- 
sure over the virtuality of the branching pseudo-parton 
ij. The parton k may or may not serve to define the 



d$(L) = d^n-^d^ij; L) , 



(2.19) 
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where ij is an on-shell parton, and 

d& 2 (i,j]Si)=d$ 2 (i,j) (2.20) 
d$ 2 (i,j;C ij )=d$ 2 (i,j) , (2.21) 

Pi\\Pj 

d^j-SC^=d^j)\^^ (2.22) 

are obtained by neglecting constant terms in the corresponding limits. Obviously, 
d<&(L) is not an exact representation of the full phase space any longer, i.e. d$(L) ^ 
<i<3>; however, this does not introduce any approximation in the procedure, since the 
subtraction counterterm is subtracted and added back in the physical cross section 
(eq. (2.2)). Finally, the asymptotic behaviours r(L) appearing in eq. (2.15) are 
directly taken from factorization formulae, eqs. (2.16) and (2.17), without any further 
manipulation. 

In the dipole formalism of ref . [9] , the following identity is exploited 

1 1 1 . 

+ 7 \ (2-23) 



Pi -PjPi-Pk Pi- PjPi ■ (Pj +Pk) Pi- PkPi ■ (Pj + Pk ) 

for the eikonal factor of eq. (2.16). The two terms on the r.h.s. of eq. (2.23) are 
symmetric for j <-> k, and thus only the first one actually needs to be considered. As 
far as collinear configurations are concerned, this term is singular only when p^ || pj, 
but not when pi || pk- Thus, the identity in eq. (2.23) has the same function as the 
partition of the phase space of ref. [8]. Furthermore, an exact parametrization of the 
phase space is chosen 

d$ = d& n - 1 fak)d& 2 {i,y,k) (2.24) 

which differs from the one of eq. (2.18) in that the parton k (called the spectator) 
plays a fundamental role, since it allows to put on shell the splitting pseudo-parton 
even if i and j are not exactly collinear, or i is not soft. Thanks to this property, in 
ref. [9] we have 

d$(L) = d$ , L — Si, C i:i , SQj, (2.25) 
and thus eq. (2.15) becomes 

fd$ = D m d$ , (2.26) 
Dij, k = r{Si, k ) + r(Cij) - r(SC tJ ) , (2.27) 

where r(Si-h) is obtained by taking the soft limit of the real matrix element, and 
keeping only the first term on the r.h.s. of eq. (2.23). 

In summary, the common feature of refs. [8, 9] is the fact that, by disentangling 
(with different techniques) the two collinear singularities that appear in each eikonal 
factor, they define subtraction formalisms based on building blocks which all have a 
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collinear-like topology; we shall denote by 1 this topology, depicted in fig. 2 (the par- 
ton labels are obviously irrelevant for topological considerations), a notation which 
is reminiscent of a branching after which the list of resolved partons is diminished 
by one unity. The explicit expressions for the building blocks, which originate from 
the soft, collinear, and soft-collinear limits and include the treatment of the phase 
space, depend on the formalism; however, in all cases they are combined according 
to eq. (2.15), i.e. according to the ^-prescription, in order to construct the sought 
subtraction counterterm. 

3. Subtraction at NNLO 

In this section, we shall introduce a general framework for the implementation of a 
subtraction method to NNLO accuracy. We shall consider the process 

e + e~ — > n jets. (3.1) 

In this way, all the intricacies are avoided due to initial-state collinear singularities, 
which allows us to simplify the notation considerably. The systematic construction 
of the subtraction counterterms that we propose in the following will however be 
valid also in the case of processes with initial-state hadrons, since the procedure 
is performed at the level of short-distance partonic cross sections. On the other 
hand, we do not present here the explicit parametrizations of the phase spaces for 
the case of initial-state partons, and we do not consider the contributions of initial- 
state collinear counterterms which are necessary in order to achieve the complete 
cancellation of infrared singularities for processes with QCD partons in the initial 
state. 

3.1 Generalities 

At NNLO, the process (3.1) receives contribution from the following partonic sub- 
processes 

e + e~ — ► m partons, m — n, n + 1, n + 2. (3.2) 
We write the amplitude corresponding to eq. (3.2) in the following way 

M m = 9 r 2 M% + g?M$ + g?+ 2 M$ + . . . , (3.3) 

where JA$, M.$ and M$ are the tree-level, one-loop and two-loop contributions 
to the process (3.2) respectively. Squaring eq. (3.3) we get: 

\M n \ 2 = gf~* \M^\ 2 + g'r 2 {M^M^ + M^M^) 

+9T (Kf + M^M^ + M 0) *M 2) ) + 0(gT + % (3.4) 
\M n+1 \ 2 = g 2 r 2 K°|i| 2 + gf {M^IMS + ^iTi-^i+i) + 0(9f + % (3.5) 



\M n+2 \ 2 = gf MZ 2 ' + 0(g 2 s n+2 ). (3.6) 



,(o) 



2 
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In eq. (3.4) the number of partons coincides with the number of jets of the physical 
process, and this implies that all partons must be resolved, i.e. hard and well sep- 
arated. On the other hand, in eqs. (3.5) and (3.6) the number of partons exceeds 
that of jets, which means that one and two partons respectively are unresolved in 
these contributions. The accuracy with which the various terms in eqs. (3.4)-(3.6) 
enter the cross section can be read from the power of g s . The Born contribution 
is proportional to gf 1 ' 4 , and appears solely in eq. (3.4). The NLO contributions 
are proportional to g 2n ~ 2 , and appear in eqs. (3.4) and (3.5). The divergences of 
the former are entirely due to the loop integration implicit in whereas those 

of the latter are obtained analytically after applying the subtraction procedure de- 
scribed in sect. 2 (there, g 2 ™~ 2 -Mn+i has been denoted by r). We understand 
that, in the actual computation of an infrared-safe observable, the matrix elements 
in eqs. (3.4)-(3.6) are multiplied by the relevant measurement functions. 

The NNLO contributions are proportional to g 2n , and we can classify them ac- 
cording to the number of unresolved partons. In the double-virtual contribution, 

2 

M$ + (M^M^* + M$*M$), all partons are resolved. The term in round 
brackets is identical to the NLO virtual contribution, except for the fact that one- 
loop results are formally replaced by two-loop ones; on the other hand, the former 
term is typologically new. However, for both the structure of the singularities is 
explicit once the loop computations are carried out. We then have the real-virtual 
contribution, Ain+i-Mn+i + -^i+i-^n+i) i n which one parton is unresolved. This is 
again formally identical to the NLO virtual contribution, but there is a substantial 
difference: in addition to the singularities resulting from loop integration, there are 
singularities due to the unresolved parton, which will appear explicitly only after 
carrying out the integration over its phase space. In order to do this analytically, a 
subtraction procedure will be necessary, and the methods of sect. 2 may be applied. 
When doing so, however, at variance with an NLO computation the analogue of the 
first term on the r.h.s. of eq. (2.2) will not be finite, because of the presence of the 
explicit divergences due to the one-loop integration. This prevents us from setting 
e — > as in eq. (2.2), and thus ultimately from computing the integral, since this 
integration can only be done numerically. It follows that the straightforward appli- 
cation of an NLO-type subtraction procedure to the real-virtual contribution alone 
would not lead to the analytical cancellation of all the divergences. 

We shall show that such a cancellation can be achieved by adding to the real- 
virtual contribution a set of suitably-defined terms obtained from the double-real 



contribution 



M-n+2 ■ This contribution is characterized by the fact that two par- 
tons are unresolved and, analogously to the case of the real contribution to an NLO 
cross section, all of the divergences are obtained upon phase space integration, a task 
which is overly complicated due to the substantial amount of overlapping among the 
various singular limits. 
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3.2 The subtracted cross section 

In this section we shall use the findings of sect. 2 as a template for the systematic 
subtraction of the phase-space singularities of the double-real contribution, which 
we shall obtain by suitably generalize the ^-prescription. In order to do this, we 
start from listing all singular limits that lead to a divergence of the double-real 
matrix elements. We observe that the S a , C a , and SC a configurations described in 
sect. 2 are also relevant to the double- real case. In addition, we have the following 
configurations [32, 33]: 

• Soft-collinear: two partons have parallel three-momenta, and a third parton 
has vanishing energy 1 . 

SC ijk : (p^O^Wph). (3.7) 

• Double soft: two partons have vanishing energy. 

SSy. (3.8) 

• Double collinear: three partons have parallel three-momenta, or two pairs of 
two partons have parallel three- momenta 2 . 

CCi jk : {pi || pj || p k ), (3.9) 
CCi jk i : (pi || pj,p k || pi). (3.10) 

• Double soft and collinear: two partons have vanishing energy, and two partons 
have parallel three-momenta. 

SSC l3 : (p?-0, P ;->0,fl||#), (3-H) 
SSC ijk : (p° - 0,p° - 0,p, || p k ). (3.12) 

• Double collinear and soft: as in the case of double collinear, but one of the 
collinear partons has also vanishing energy. 

SCCi jk : (p°i -> 0,pi || pj || p k ), (3.13) 
SCC ijH : (p°i -> 0,Pi || Pj,Pk || Pi)- (3.14) 

1 Note therefore that, upon removing the parton labels, the SC symbol denotes both the cases 
with one or two unresolved partons, eqs. (2.6) and (3.7) respectively. 

2 The former case is usually denoted as triple collinear. We prefer this notation since, if a \\ b 
and b \\ c, necessarily a \\ c. Furthermore, the present notation is more consistent with the strongly- 
ordered limit case. It gives minimal but sufficient information. 
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• Double soft and double collinear: as in the case of double collinear, but two of 
the collinear partons have also vanishing energy. 



The soft and collinear limits in each of eqs. (3.7)-(3.16) are understood to be taken 
simultaneously, following for example the rules given in ref. [33] (see in particular 
eqs. (23) and (98) there). 




Figure 3: NNLO-type topologies: 2 (left panel) and 1<8> 1 (right panel). 

The necessity of introducing the notion of topology emerges at NNLO even with- 
out considering the soft limits and the problem of overlapping divergences, as is clear 
by inspection of the purely collinear limits, eqs. (3.9) and (3.10). These are asso- 
ciated with the branching processes depicted in fig. 3; we denote the corresponding 
topologies by 2 and 1(8)1 respectively (again, this notation serves as a reminder of 
the number of partons to be removed from the list of resolved partons). Some of 
the singular limits in which one or two partons are soft cannot be straightforwardly 
associated with either topology. However, as in the case of NLO computations, this 
can be done after some formal manipulations, whose nature (be either a partition 
of the phase space, or a partial fractioning, or something else) we do not need to 
specify at this stage. Suffice here to say that, after such manipulations, all singu- 
lar limits in eqs. (3.7)-(3.16) that feature at least one soft parton will in general 
contribute to both 2 and 1 <g> 1 topologies (we may formally write, for example, 
f(SS) = f(SS^) + f(SS^)). The same kind of procedure can be applied to 
the singular limits of eqs. (2.4)-(2.6), since topology 1 can always be seen as a sub- 
topology of 2 or 1(81; thus, S, C, and SC limits will be manipulated, if need be, 
so as to be associated with topologies 2 and 1<8>1. In the case we shall need to dis- 
tinguish between the singular limits in the various topologies, we shall denote them 
by L( T \ with T = 1, 2, 1(8)1. However, as for parton labels, topology labels may be 
understood in what follows. 

We now claim that by applying the ^-prescription defined in sect. 2 we can 
systematically subtract the singularities of the double-real contribution to the NNLO 



SSCCijk : 
SSCCijki : 



(Pi -* °>P? -»■ °>ft II Pj II Pk), 
(Pi -»■ 0>P° -* °>Pi II Pj>Pk II Pi)- 



(3.15) 
(3.16) 
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cross section, provided that the operator S NLO is replaced by £ NNLO , where 

s NNLO (f) = s NLO (f) - E E {/ (sc<p) + / (ssp) + / {ccp) 

+ / (SSC<P) + / (SCCP) + / (SSCCP) j, (3.17) 

and a denotes all indices relevant to the corresponding singular limits, which can 
be read in eqs. (3.9)-(3.16). As in the case of NLO computations, the iterative 
^-prescription comes to an end thanks to the infrared safety of the observables. 
However, at NNLO up to four iterations are necessary in order to define /, which 
can be easily generated by means of an algebraic-manipulation code. In such a way, 
one obtains up to 51 terms for each topology; fortunately, such a massive counterterm 
can be greatly simplified. We start by observing that the freedom in the definition 
of the asymptotic form of the matrix elements associated with a given singular limit 
allows us to exploit commutation properties as done in eq. (2.13). Furthermore, the 
presence of S NLO in the definition of S NNLO implies that some of the terms obtained 
with the ^-prescription will be formally identical to those appearing in an NLO 
subtraction. This suggests us to write 

f = f-i + fS + fT\ (3-18) 

where 

/-i = E / (sf } ) + E / - E / , (3.19) 

i i<j ij 

and 

= E {/ (ccf>) + / (ssp) - f (sc<p) -f([c® ss\p) 

a { 

- f ([CC © s\<p) - f ([CC © SS}^) -f([S® s\<p) - f ([C © C}^) 
+ f([C®S® S}^) + / ([CC © s © s\<p) + / ([C © C © S}^) 

+ f([C®C® SS\<P) -f([C®C®S® 5]f >) j, (3.20) 

where we used the fact that (for example) SSC = SS © C = C © SS, and a denotes 
symbolically the relevant parton indices, not indicated explicitly in order to simplify 
the notation. The physical meaning of eqs. (3.18)-(3.20) is clear: in the double real 
contribution to an NNLO cross section, there are terms with one or two unresolved 
partons. The former have the same kinematics as those relevant to a pure NLO 
subtraction. More interestingly, they also have the same kinematics as the real- 
virtual contribution. Although this fact formally results from the application of the 
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^-prescription, it can also be understood intuitively: by requiring more stringent jet- 
finding conditions (for example, by enlarging the minimum p T which defines a tagged 
jet), the NNLO cross section turns into an NLO one, which receives contributions 
only from the real-virtual term, and from the pieces obtained by applying S NLO to 
the double-real matrix elements. Clearly, if / is associated with an (n + 2)-body 
final state, each term in /_! and /_ 2 factorize (n + l)-body and n-body hard matrix 
elements and measurement functions respectively. 

In order to give more details on the structure of the subtraction that emerges 
from the ^-prescription, let us denote by 

w = g? (K 1} f + M^M^ + M^M^) , (3.21) 
rv = gf (M$.MZ + M^Mnl) , (3.22) 

(3.23) 



2n 

rr = g s 



M 



(0) 
n+2 



2 



the double- virtual, real- virtual, and double-real matrix elements squared respectively, 
possibly times measurement functions that we understand. The jet cross section is 

da = da rr + da rv + da vv (3.24) 
= J rrd$ n+2 + f rvd$ n+1 + J vvd$ n . (3.25) 

We start by applying the ^-prescription to the double-real contribution 

da rr = J (rrd<S> n+2 - m^ 2 d^ 2 +2 - rr_id$-| 2 )+ J rr_ 2 d^~l 2 + J rr^d®'^ (3.26) 

where we allowed the possibility of adopting two different parametrizations for the 
phase spaces attached to rr_ 2 and to fr_ 1 . The first integral on the r.h.s. of eq. (3.26) 
is finite, and can be computed numerically after removing the regularization by 
letting e — > 0. The second and the third integrals will contain all the divergences of 
the double-real contribution, to be cancelled by those of the real-virtual and double- 
virtual contributions. However, only the second term can, at this stage, be integrated 
analytically. In fact, ff_i is, according to its definition, obtained by considering 
the asymptotic behaviour of the double-real matrix element squared in the singular 
limits S, C, and SC of NLO nature. These limits will render manifest only part 
of the singular structure of the matrix elements, preventing a complete analytical 
integration of the divergent terms. For example, if parton % becomes soft, rr(Si) 
factors the eikonal term of eq. (2.16), and the integration over the variables of parton 
% can be carried out analytically, as in the case of NLO. However, at NNLO we must 
take into account that another singularity may appear - say, parton I may also 
become soft, and rr(Si) has too complicated a dependence upon the variables of 
parton I to perform the necessary analytical integration. Notice that this is not true 
for, say, rr(SSu), which is why the second term on the r.h.s. of eq. (3.26) can indeed 
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be integrated. This suggests combining rr_i with the real- virtual contribution which, 
as discussed in sect. 3.1, cannot be integrated analytically too (although for different 
reasons). In order to achieve this combination, we must choose the phase-space 
measure appropriately; in particular, we shall use 

d$£ 2 = d$ n+1 d^ 1 , (3.27) 

where d<f> n+ i is the exact (n + l)-body massless phase space, and dQ 2 l is related to 
the phase space relevant to the partons whose contributions to the singular behaviour 
of the double- real matrix elements is explicit in ff_i, and includes the measure over 
the virtuality of the branching parton; we shall show in what follows how to construct 
explicitly the phase spaces of eq. (3.27). Having done that, we define a subtracted 
real-virtual contribution as follows 

rv {s) =rv + J ff_id$2 1 ■ ( 3 - 28 ) 

In an NLO computation, eq. (3.28) would amount to the full NLO correction, with 
the first and the second term on the r.h.s. playing the roles of virtual and real 
contributions respectively. This implies that the explicit poles in 1/e that appear in rv 
because of the loop integration will be exactly cancelled by those of rf-i which result 
from the phase-space integration d^J 1 . Thus, although in an NNLO computation 
rv ^ still contains phase-space divergences, we can manipulate it in the same manner 
as a real contribution to an NLO cross section: 

J rv^d$ n+1 = J (™«d$ n+1 - H; (s) rf$ n+ i) + J ™ (s) di» n+ i . (3.29) 

The last term on the r.h.s. of eq. (3.29) can now be integrated analytically, whereas 
the first is finite and can be integrated numerically. Therefore, by combining part 
of the double-real and the real-virtual contributions we have managed to define a 
scheme in which the analytical integration of all the divergent terms is possible. By 
combining eqs. (3.25), (3.26), and (3.29) we get 

da = J (rrd<$> n+2 - ff„ 2 rf$^ 2 - rr^d^'^) 
+ J (rv«d$ n+1 - ™«d$ n+1 ) 

+ Jrf- 2 d^~l 2 + j rl) {s) d$ n+1 + Jvvd$ n , (3.30) 

where the first two terms on the r.h.s. are finite, and can be integrated numerically 
after letting e — * 0; the sum of the remaining terms is also finite, but they are 
individually divergent and must be computed analytically. 

In order to show explicitly how our master subtraction formula eq. (3.30) works, 
we consider the unphysical case in which only singular limits of collinear nature 
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can contribute to the cross section. In such a situation, the application of the E- 
prescription is trivial, and we readily arrive at 3 : 

rJ-_ 2 d^~l 2 = rr{CC) d^~l 2 {CC) - rr{C © C) d®- 2 +2 {C © C) , (3.31) 
rr- 1 dQ£ 2 = rr(C)dQ£ 2 (C). (3.32) 

As in the general formula for the subtraction counterterm at NLO, eq. (2.15), we 
leave the possibility open of associating different phase spaces measures with different 
terms in fr_ 2 . This has the advantage that each parametrization can be tailored in 
order to simplify as much as possible the analytical integration. The drawback is 
that, in general, there could be algebraic simplifications among the contributions to 
rr_ 2 (here, between rr{CC) and rr{C ©C)), which can be explicitly carried out only 
upon factorizing the phase space. We shall discuss this point further in the following, 
in the context of a more physical case. 

By construction, rr(C © C) is the strongly-ordered, double-collinear limit. This 
may coincide with rr(CC), and it does in particular in topology 1<8>1, but in general 
fr_ 2 is different from zero 4 , and corresponds to the non-strongly-ordered part of 
the double-collinear limit. As such, when a kinematic configuration is generated in 
which two, and only two partons are collinear, we have fr_2 — > 0; in other words, 
the C limit of fr_ 2 is zero. This is what should happen: in the collinear limit, rr(C) 
(which appears in fr_i) is sufficient to cancel locally the divergences of rr, and thus 
fr_2 should not diverge in this limit, since otherwise the first term on the r.h.s. of 
eq. (3.30) would not be finite. Analogously, in the double-collinear limit the local 
counterterm for rr is rr(CC); therefore, in order to avoid divergences, in such a limit 
the contribution of rr{C © C) (in fr_ 2 ) must cancel that of rr(C) (in rr_i). Notice 
that, for these cancellations to happen not only at the level of matrix elements, but 
also at the level of cross sections, suitable choices of the phase spaces must be made. 
The C limit thus relates d&~+ 2 (CC) to d&~+ 2 (C © C), whereas the CC limit relates 
d$~+ 2 (C ©C) to d<&~l 2 (C). A suitable choice in the former case clearly includes the 
trivial one, where the two phase spaces are taken to be identical. 

We finally note that the subtraction formula of eq. (3.30) appears in a very 
similar form in ref. [21], where a first discussion was given on the extension of the 
subtraction method to NNLO. Ref. [21] constructs the NNLO subtraction formula 
building upon the NLO dipole subtraction [9], and provides the explicit expressions 
of the NNLO counterterms for the leading-colour contribution to e + e~ — > 2 jets. 
Ref. [21], however, does not discuss the way in which the subtraction kernels can be 
constructed for more general processes. In our approach, the ^-prescription provides 
a general framework for the construction of the counterterms rr_i, rr„ 2 and rv^ s \ 
and explicitly suggests the subtractions of eq. (3.30). We find it reassuring that we 

3 The partem indices play an obvious role here, and we omit them in order to simplify the notation. 
4 Hcre and in what follows, "zero" means non-divergent. 
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arrive at a subtraction structure consistent with that of ref. [21], given the fact that 
neither here nor in ref. [21] a formal proof is given that eq. (3.30) achieves a complete 
cancellation of the infrared singularities. We do obtain such a cancellation explicitly, 
for the C F T R colour factor of the dijet cross section, as we shall show in the next 
section. On the other hand, to the best of our knowledge the counterterms presented 
in ref. [21] have not yet been integrated over the corresponding phase spaces. 

4. An application: the C F T R part of e + e — > 2 jets 

In this section, we shall apply the subtraction procedure discussed in sect. 3 to a 
physical case, namely the contribution proportional to the colour factor C F T R of the 
dijet cross section in e + e~ collisions. This is a relatively simple part of the com- 
plete calculation of an observable to NNLO, but it allows us to discuss the practical 
implementation of most of the features of the subtraction procedure we propose in 
this paper. The branching kernels that we shall introduce are universal, i.e. they 
can be used in any other computations where they are relevant. We shall also define 
precisely the phase spaces needed to integrate the above kernels over the variables 
of unresolved partons, i.e. the quantities d<& used in the formal manipulations of 
sect. 3. We shall show that the subtraction procedure leads to the expected KLN 
cancellation, and that the numerical integration of the finite remainder gives a result 
in excellent agreement with that obtained in ref. [19]. 

Although not necessary for the implementation of eq. (3.30), in this section 
we shall use partial fractioning to deal with the eikonal factors associated with soft 
singularities, and combine them with the corresponding collinear factors by using the 
same kinematics in the hard matrix elements that factorize. Since such kinematics 
will also enter the measurement functions appearing in the counterterms, all of the 
manipulations involving the hard matrix elements will also apply to the measurement 
functions; for this reason, the latter will be left implicit in the notation. 

We shall reinstate here the notation commonly used for QCD amplitudes, which 
can be written as vectors in a colour space including the coupling constant. Thus, 
eq. (3.3) is now rewritten as 

\M aia2 ... am ) = |-Mi!... a J + |-M£l 2 ...aJ + lM!...aJ + • • • (4-1) 

where a\a2 . . . a m are flavour indices, and particles other than QCD partons are 
always understood. When not necessary, flavour labels and colour vector symbols 
may also be understood. 

UV renormalization is performed in the MS scheme, just by expressing the bare 
coupling a" in terms of the renormalized coupling a s (fj, 2 ) at the renormalization scale 
\i. We use the following expression 

a u s nl e S e = a s (fi 2 )fi 



_ ajji ) f3o 
2tt e 



(4.2) 
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where (3q = ^-C A — ^n F T R , and 

S e = (47r) e e- £7B (4.3) 

is the typical phase space factor in d = 4 — 2e dimensions (7^ = 0.5772... being the 
Euler number). 




Figure 4: Sample diagrams for double-real (left panel) and real-virtual (right panel) 
contributions to the C F T R part of the dijet cross section. 



The diagrams contributing to the C F T R part of M.f^ (i.e., to the double real) are 
of the kind of that displayed in the left panel of fig. 4; thus, the singular limits we are 
interested in are associated with the branchings g — > qq and q — > qq'q' (the identical- 
flavour branching q — > qqq is also accounted for by considering q — > qq'q', since its 
interference contributions are not proportional to C F T R ). As far as the real- virtual 
contribution is concerned, a sample diagram of A4^ is depicted in the right panel 
of fig. 4; the singular limits of phase-space origin which, according to the discussion 
given in sect. 3.2, are obtained upon applying the ^-prescription, are associated with 
the branching q^> qg. 

4.1 The double-real contribution 

Thanks to the universality properties of soft and collinear emissions, any matrix 
element with singularities due to the branchings g — > qq and q — > qq'q' can be used 
to define the subtraction terms relevant to the double-real contribution. Thus, we 
write the double-real matrix element squared as 

rr = {M%, q ^ an JM^ q ^ an J , (4.4) 

where the labels imply that the four- momenta of q', q' and q are ki, A; 2 and k% 
respectively. When applying the ^-prescription to eq. (4.4) we obtain from eqs. (3.19) 
and (3.20) 

ff_ 2 = rr{CC) + rr(SS) - rr(CC © SS) 

- rr{C © C) - rr(C © SS) + rr{C ®C®SS), (4.5) 
fr_i = rr(C) , (4.6) 
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where clearly C = C\ 2 , CC = CC123, and so forth. In order to compute explicitly the 
quantities that appear in eqs. (4.5) and (4.6), we shall use the results of ref. [33] (see 
also refs. [32, 34]). The asymptotic behaviours obtained through successive iterations 
of S NNLO , characterized by the © symbol, can be freely defined to a certain extent. 
We give such definitions within a given parametrization of the phase space, which 
we shall introduce in the next subsection. 



4.1.1 Choices of phase spaces 

As discussed in sect. 3.2 (see in particular eq. (3.26)), the definition of the double- 
real subtraction terms implies the necessity of defining d&~+ 2 and which are 
related to some extent. As a preliminary step, we make here the choice of associating 
the same d<&~+ 2 with all of the terms that appear in eq. (4.5); this is in principle 
not necessary (see eq. (3.31)), but we find it non restrictive in the computations that 
follow (more complicated kernels may require a different choice). We then write the 
exact (n + 2)-body phase space using eq. (B.ll) 

d<S> n+2 = ci$ n (123) d$ 3 (l, 2, 3) , (4.7) 
Zir 

where we use the momenta labels rather than the four-momenta to shorten the 
notation (consistently with eq. (2.18), 123 means k\ + k 2 + fc 3 ), and 

S123 = (h + h + hf . (4.8) 

As in the case of NLO computations, eq. (4.7) is unsuited for the analytical integra- 
tion over the variables of unresolved partons since 123, which enters the phase space 
associated with the non-singular part of the matrix element, has an off-shellness S123, 
while the reduced matrix element which corresponds to such non-singular part has 
all of the final-state QCD partons with zero mass. We follow the same strategy as 
in eq. (2.19), i.e. that of ref. [8]: we introduce a four- momentum 123 with invariant 
mass equal to zero 5 , and define 

d§- 2 +2 = rf$ n (123) ^ d$ 3 (l, 2, 3) . (4.9) 

This definition is equivalent to that in eq. (2.19), and differs from that in eq. (2.24) 
in that no spectator is used to keep the branching parton on shell. 
For the definition of d^~]_ 2) we use again eq. (B.ll) 

d® n+ 2 = ^ d^ n+1 (12) d$ 2 (l, 2) , (4.10) 



5 123 can be denned to have the same three-momentum as 123 in the e+e rest frame, and zero 
mass. However, the precise definition is irrelevant in what follows. 
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from which the analogue of eq. (4.9) can readily follow: 

dsr. 
2tt 



d$£ 2 = d$ n+ i(12) ^d$ 2 (l, 2) . (4.11) 



There is however a subtlety: as discussed in sect. 3.2, the choice of d<b n \ 2 must match 
that of d&~ 2 2 . m order to see how this can happen, we use eq. (B.22) to write 

d$£ 2 = d$ B (123) ^ ^ d<t> 2 (3, 12) d$ 2 (l, 2) . (4.12) 

The part relevant to the 12 — > 1 + 2 branching is identical in eqs. (4.11) and (4.12), 
except for the fact that the invariant mass of the 12 system, s i2 , enters <i<3> 2 (3, 12), 
and not only d$ 2 (l,2), in eq. (4.12). If this dependence could be safely neglected 
in the limit s i2 — > 0, i.e. in the singular region in which d&~]_ 2 an d d&~ +2 must 
match, the choice of eq. (4.11) would be appropriate. However, as can be seen from 
eqs. (B.25) and (B.27), rf$ 2 (3, 12) contains the factor 

(z 3 z 12 s 123 - z 3 s 12 )~ e , (4.13) 

which acts as a regulator in the integrations in <i<I> 2 (3, 12); since these are in general 
divergent, the regulator affects finite (and possibly also divergent) terms, and thus 
cannot be ignored even if the limit s i2 — > is considered. The regulator of eq. (4.13) 
is implicit in d<& n+ \(12) in eq. (4.10); the effect of neglecting it as done by defining 
d$~} 2 in eq. (4.10) therefore leads to neglecting contributions to the cross section if 
d$ n+ i(12) is involved in the integration of divergent terms. This is what happens, 
since the subtracted real- virtual contribution is indeed divergent. This problem may 
seem to be cured by inserting the regulator of eq. (4.13) into the r.h.s. of eq. (4.11). 
However, the definition of g?$~+ 2 should not depend on whether the branching 12 — > 
1 + 2 is followed by the branching 123 — > 3 + 12, which is what the presence of 
Si 23 in eq. (4.13) implies. The most general form for the regulator to be inserted in 
eq. (4.11) can be deduced by writing 

d$ n+2 = ^ dg3 : (n+2) d$ 2 (12, 3 . . . (n + 2)) d$ 2 (l, 2) d$ n (3, . . . , (n + 2)) . (4.14) 
The first phase space on the r.h.s. features the regulator 

{z\2Z 3 ...( n +2)Q 2 - Z 12 S3...( n +2) ~ ^3...(n+2) S12) = (^3...(n+2)sST) ^1 ~ ) 

(4.15) 

which can be obtained from eqs. (B.17) and (B.18), and where 

S%£* = (zi2Z3...( n +2)Q 2 - Z 12 S 3 ... (n+ 2)) /^3...(n+2) , (4-16) 

with Q 2 = (ki + . . . + k n+2 ) 2 . Eq. (4.15) suggests to replace eq. (4.11) with 

d4&2 = ^ n+1 (12) ^ (l - ^) d* 2 (l, 2) , (4.17) 
where the quantity s™ 2 x is a constant with respect to the integration in s 12 . 
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4.1.2 Computation of the divergent terms 

In this section, we integrate rr_ 2 and rr_ 1 given in eqs. (4.5) and (4.6) over the phase 
spaces defined in sect. 4.1.1. We start from the case of ff-2, and write 



rr_ 2 



q...a n+ 2 



(4.18) 



(3) 

The kernel K ql -, q is a matrix in the colour space, and collects all singular behaviours 
associated with the branching q — > q[q' 2 q 3 according to the combination given in 
eq. (4.5). For its explicit construction we shall use the factorization formulae given 
in ref. [33]. In general, we shall denote the branching kernels as follows 



K (n) 



a\...a n 



for S(ai, . . . , a n ) — > cti + . . . + a r 



(4.19) 



where S(ai, . . . , a n ) is the parton flavour which matches the combination a± + . . .+a n . 
In this paper, only n = 2, 3 will be considered. Notice that the notation of eq. (4.19) 
does not coincide with the one typically used for two-parton branchings, where the 
Altarelli-Parisi kernel is denoted by P^si^^)- There is a remarkable simplification 
that occurs in the computation of eq. (4.18); namely, one can prove that (see app. A) 



rr(SS) = rr(C ® SS) , rr(CC ® SS) = rr(C ®C®SS). 



(4.20) 



Thus 



7f_ 2 =rr{CC)-rr{C®C). 



(4.21) 



Since in the collinear limits colour correlations do not appear, eq. (4.21) implies that 

(3) 

K,,, F ,,„ has a trivial structure in colour space. Using the same normalization in the 

11I2I3 

factorization formulae as in ref. [33] (see eq. (29) there), we thus rewrite the kernel 
as follows 

,s "'"""' l2 KS,-^S,). <«■*) 



6 123 



where the first term can be read from eq. (57) of ref. [33] 

4z 3 + ( Zl - z 2 f 



K cc 



and 



1 n rp S 123 



t 2 

t 12,3 



Sl2Sl23 



+ 



Zl + Z 2 



+ (1 - 2e) ^1 + z 2 - 



S12 



s 123 



ZiSjk ZjSik 
Z{ -\- Zj 



Zl "T Zj 



it] ■ 



(4.23) 
(4.24) 



The quantity in eq. (4.22) corresponds to the successive branchings q — >• q^g, 

g — > q[q' 2 . We can therefore construct it by using the polarized Altarelli-Parisi kernels 



Pjfq and , as shown in app. A. The result is 



K C®C _ 1 r rp «123 

^a'a'a ~^ Ly F 1 R 



A% | 4z 3 + ( Zl - z 2 f | 
S12S123 z x + z 2 



2e) ( Zl + z 2 ) 



(4.25) 
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where 



(5) = 2 gigjfc (4 ^ 2g) 



ij ' fc ^ + Zj 



We note that the result in eq. (4.25) corresponds to the naive s± 2 — > limit of K^,^ 
in eq. (4.23), since £12,3 — > t[ 2 \ oc v /si2 in such a limit (see app. B). The kernel 

'(3) ' 

the results of eqs. (4.23) and (4.25), we find 



q'^Qs * s therefore manifestly simpler than its two contributions in eq. (4.22); using 



(3 ) _ (87ra>g e ) 2 1 s 123 
**** ~ sj 23 2 Cf1r — 2 



t 2 - t iS)2 

'-12,3 H2,3 _ (I _ 2c) 2 



(4.27) 



512^123 s 123_ 

We can now obtain the analytical expressions of the 1/e poles in ff_ 2 . Using eq. (4.9) 

/ rr. 2 d*£ 2 = \M q °l n+2 \ 2 d$ n J ^^^$3(1,2,3) . (4.28) 

We use the parametrization of eq. (B.33) to perform the integration in eq. (4.28). The 
invariants are expressed in terms of the integration variables according to eqs. (B.30) 
and (B.31); however, one must be careful when replacing these expressions into 
since in doing so finite terms are generated when s 12 — > 0; such terms should not 



appear, since they have been explicitly neglected when working out. To avoid 

this, we use 



Z1S23 - z 2 s 13 -> v/s^o , 6 = 2^3(2(1 - (2)^123(1 - z 3 )x (4.29) 
(s) 

in the computation of t\ 2 3 - In a more general case, the replacements 

S13 -> Si23(l - C2) , S23 = S123 C2 (4.30) 
should also be made. The result of the integral in eq. (4.28) is 

/ tXk, «.(!• 3) = m 2 C F T» (m 2 ' i (! + f «) + 0(£ , , (4.3!) 



where we denoted by s" 2 | x the upper limit of the integration in S123, whose form does 
not need be specified here; a definite choice will be made in sect. 4.3. 
We now turn to the case of rf_i, and write the analogue of eq. (4.18) 

rr-, = (M^ 2 \K^\M^ a ) • (4.32) 

As in the case of K^-,^, the kernel -K"^ is purely collinear (see eq. (4.6)), and 
therefore has a trivial colour structure. Using again the factorization formulae with 
the normalization of ref. [33] (see eq. (7) there), we have 

M<f 5 .M = ^^r, (4-33) 
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where P q q is given in eq. (A. 8). Spin correlations are essential for the kernel to 
be a local subtraction counterterm, but they do not contribute to the analytic inte- 
gration which follows, where the spin-dependent collinear kernel P q q can be replaced 
by 



1 - 



2z(l-z) 
1 -e 



(4.34) 



Using eq. (4.17) we obtain 



/ 



M (0) 

3<j3...a n+ 2 



The integral can be performed straightforwardly, with the result 
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2tt W> 
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max 
12 , 
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3 + Y e 



'56 _ tt2 
27 9 
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81 27 n ~ " 3 Cs ] 63 



x 1-^-2^3 +0(e 3 ), 



(4.36) 



where all the terms up to C(e 2 ) have been kept, since eq. (4.36) will be used to define 
the subtracted real- virtual contribution as given in eq. (3.28), which will generate 
further poles up to 1/e 2 . 



4.2 The real-virtual contribution 

In this section, we shall construct the subtracted real-virtual contribution, defined 
in eq. (3.28). As discussed there, rv^ has no explicit 1/e poles, since those resulting 
from the one-loop integrals contributing to the (unsubtracted) real-virtual contribu- 
tion rv are cancelled by those of rf_i. We can show this explicitly for the C F T R part 
of the e + e~ dijet cross section: the contribution of 

rv = (Mf q \\M^ q \) + (M$ q \Mf q \) (4.37) 

to this colour factor is entirely due to the T R part resulting from UV renormalization: 

a s 2T R n F (0) 2 

TV \t r = ^3^— M 9Q9 ■ ( 4 - 38 ) 

By using eqs. (3.27), (4.35), and (4.36), with gq 3 . . . a n+2 = gqq and summing over 
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n F quark flavours, it is apparent that rv^ is indeed free of explicit 1/e poles: 



rv 



rv\ TR d$ 3 + J rr^d^ 4 l 

flav 

f 2 ( /i 2 V 



a s T R n I 
2ir e 



M 



(o) 
ggg 



„max 
-12 , 



r(i-c) 



2 10 /56 7T 

3 + y 6+ 27 "T"'" 



/328 5 2 4 \ o 
— 5 + 3l0g ' — 



„max 
6 12 , 



7T 



1 - -e" - 2e% ) \ + 0(e 3 ) 




M 



(o) 
ggg 



d§ 3 + C(e) . 



(4.39) 
(4.40) 



We underline the presence of s™* in eqs. (4.39) and (4.40), due to the integration 
performed in eqs. (4.35) and (4.36), and to the choice of phase space of eq. (4.17). Its 
specific form depends on the branchings involved in the construction of fr_i, and can 
be read from eq. (4.16); we shall write it explicitly in the following, after choosing a 
parametrization for the phase spaces c?$„ +1 and d<t n+1 . 

Implicit divergences remain in eq. (4.39), whose analytic computation requires 
the definition of the subtraction counterterm rv^ which appears in eq. (3.29). It is 
clear that these divergences are those of the matrix element M g % which factorizes 



in eq. (4.39) and thus are due to the branching q 
full generality, we consider 



qg. In order to study them with 



rv 



(o) = ^(o) 



gq...a n+ i 



(4.41) 



from which we construct the subtraction counterterm by applying the ^-prescription: 



rv 



(0) 



rv^(C)+rv^(S)-rv^(SC) 



,(o) 



(0), 



(4.42) 



Clearly, rv^ is obtained by multiplying rv (0) by the r.h.s. of eq. (4.39), divided by 

M g % d&3, the same manipulations can be carried out on the r.h.s. of eq. (4.40), 
in this way obtaining the local subtraction counterterm to be used in the numerical 
computation of the integral that appears in the first term on the r.h.s. of eq. (3.29), 
or in the second line on the r.h.s of eq. (3.30). We stress that eq. (4.42) coincides with 
eq. (2.14), since the subtraction of the phase-space singularities of rv ^ is identical to 
that performed in the context of an NLO computation. In what follows, similarly to 
what done in sect. 4.1, we shall choose the same parametrization of the phase space 
for the three terms on the r.h.s. of eq. (4.42), whose sum will therefore be essentially 
identical to one of the dipole kernels of ref. [9]. As before, we shall however not use 
the dipole parametrization for the phase space, and our integrated kernel will thus 
be different from that of ref. [91. We obtain 



rv 



(0) 



V l q...a n+ i \ qg r l q...a n+1 , 



(4.43) 
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where the kernel has now a non-trivial colour structure 



K$=- Sj ^Y.<*T,-T t . (4.44) 

Here 



S 19 k+q 



V® = , 2fVP * + (1 - z)(l - e), (4.45) 
Pq ' Pk (4.46) 



(Pq + Pg) -Pk ' 

and in order to obtain the result in eq. (4.45) we have decomposed the eikonal factor 
that appears in rv^°\S) as shown in eq. (2.23) 6 . 

According to our master subtraction formula, eq. (3.30), we have now to integrate 
rv^ over the phase space, and thus we have to choose a parametrization for d& n+ \. 
A possibility is that of using again the form of eq. (4.17), which with the parton 
labeling of the case at hand would read as follows 

d$ n+1 = d* n {qg) ^ (l - ^) **2(?, 9) ■ (4-47) 

On the other hand, the present case is simpler than that discussed in sect. 4.1. The 
integration of rv^ over d$> n+ i will not be followed by another analytical integration; 
thus, the regulator that appears explicitly in eq. (4.47) is not actually necessary. Fur- 
thermore, its presence would require the presence of an analogous regulator for S123 in 
eq. (4.9). Thus, in order to simplify as much as possible the analytic computations, 
we rather adopt the following form 

ds 

d$ n+l = d$ n (qg) -7T- d ^ 9) ■ (4.48) 
Having made this choice, we can deduce the expression of s^ ax from eq. (4.16): 

*Sr = (!-*)*«■ (4-49) 

In order to integrate the kernel K^J over ds qg d§2(q, g), we use eq. (B.17), choosing 
the reference four- vector n there to coincide with pt; in such a way, the variable z\ 
of eq. (B.14) coincides with z of eq. (4.46). We also define 

Q 2 q9 k = (P q +P 9 +Pk) 2 , (4.50) 



6 We note that in this calculation the kernel needed to construct ru^ coincides with the NLO 
one owing to eq. (4.38). In more general cases, a formula similar to (4.42) for fv^ still holds, but 
rv^(C), rv (s \S) and rv^ s \SC) will have to be computed using the results of refs. [36, 37, 38]. 
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which is treated as a constant during the integration. We obtain 



r(i 



E 



qgk , 



qg,k 



In eq. (4.39) we also need 
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(0) / A* 
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*12 , 



x (MW\T q .T k \M®)d* 



(4.51) 
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\T q -T k \M q l n+1 )d<S> n , 



with s™, 3 "* given in eq. (4.49), and 
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,Q 2 
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top* f^d* 2 (q,g)V, 

J Sqg 



I% = *Kff^d* 2 (<l,9) (£) v, 

J 6 <?9 \ 6 12 / 



(2) 

2 



r(i-e) VQ 



(2) 



(4.52) 

(4.53) 
(4.54) 



A somewhat lengthy computation returns the following results 
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(4.55) 
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+ log 2 (y max ) ( - - 4 log(l - y max ) + - log(j/ n 
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+ 0(e 2 



(4.56) 



where we have introduced the arbitrary parameter y max , which defines the upper limit 
of the s qg integration 

s qg < y ma x Q\g k , < y max < 1 . (4.57) 
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The physical results must not depend on y max , whose variation amounts to changing 
finite contributions to the subtraction counterterms; whether this condition is fulfilled 
represents a powerful check on the correctness of the subtraction procedure. The 
parameter y max is the analogue of the free parameters ^ cut , 5 r and 5 Q of ref. [8]. 

4.3 Results 

In this section, we shall use the results obtained in sects. 4.1 and 4.2 in order to 
implement our master subtraction formula eq. (3.30). As a preliminary step, we 
need the C F T R part of the double-virtual contribution 



vv 



(M%\M { 3) + (M%\M*S) 



(4.58) 



which can be straightforwardly obtained from ref. [39] : 
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(4.59) 



where Q 2 is the e + e~ cm. energy squared. The subtraction counterterms in the first 
and second integrals on the r.h.s. of eq. (3.30) are constructed using the kernels 
and defined in the previous sections. It must be stressed that we shall have to 
subtract twice, since we have two configurations contributing to topology 2: the 
first is that depicted on the left panel of fig. 4, the second being identical except for 
the fact that the gluon is attached to the other quark leg emerging from the 7 — > qq 
branching. This also implies that the double-real contribution to the last line on the 
r.h.s. of eq. (3.30) will be obtained upon multiplying by two the result presented 
in eq. (4.31) (an overall factor of n F appears when summing over flavours). Finally, 
the largest kinematically-allowed value of the virtuality of the branching parton in a 
1 — > 3 splitting is Q 2 . For consistency with what done in eq. (4.57), we thus set 



max _ q1 
a 123 — i/max V 



(4.60) 



as the upper limit of the integration over S123 performed in sect. 4.1.2. The real- 
virtual contribution / rv^d$ n+ i can be obtained from eqs. (4.39), (4.51) and (4.52); 
for the dijet cross section, q . . . a n+1 = qq, and the colour algebra is trivial: 



(M 



qq I -L q 



(0)\ 



■T q \M qq 



-c F 



M 



(4.61) 



In the computation of the real-virtual contribution, we also necessarily have Q 2 qg k = 
Q 2 (see eq. (4.50)). When putting this all together, we obtain what follows for the 
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analytic-computed part of the subtraction formula of eq. (3.30): 
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(4.62) 



This result is non-divergent, which proves the successful cancellation, as dictated by 
the KLN theorem, of the soft and collinear divergences which arose in the interme- 
diate steps of the computation. 

The cancellation of the divergences does not guarantee that the finite parts result- 
ing from the implementation of the subtraction procedure are all included correctly. 
In order to check this, we started by computing the C F T R part of the NNLO con- 
tribution to the total hadronic cross section, and compared it with the well-known 
analytic result, which can be read from [11] 
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where 



R = 



a e T e" 



+ 40 



hadrons) 



+ 0{al 



(4.63) 



(4.64) 



cr(e + e — > ) 

Although the physical result must not depend on y max , the three lines on the r.h.s. of 
eq. (3.30) (which we call 4-parton, 3-parton, and analytic contributions respectively) 
separately do. We thus compute the rate for different choices of y ma , x . The results 
are presented in table 1. As can be seen there, the dependence on y max of the three 
contributions to our subtraction formula is fairly large. However, such a dependence 
cancels in the sum, within the statistical accuracy of the numerical computation. We 
defined 

I total — exact I 



pull 

exact = 



error 



11 



+ 4C 3 = -0.69177, 



(4.65) 
(4.66) 



the error being that due to the numerical computation, reported under the column 
"total" in table 1. We note that the 4-parton result changes sign when y max is 
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2/max 


4-parton 


3-parton 


analytic 


total 


pull 


0.2 


-0.692 ± 0.002 


1.4241 ± 0.0004 


-1.4236 


-0.692 ± 0.002 


0.11 


0.4 


-0.229 ± 0.004 


9.8001 ± 0.0003 


-10.2610 


-0.690 ± 0.004 


0.44 


0.6 


0.074 ± 0.02 


12.3440 ± 0.0003 


-13.0753 


-0.657 ± 0.02 


1.74 


0.8 


0.235 ± 0.005 


13.3759 ± 0.0003 


-14.2989 


-0.688 ± 0.005 


0.75 


1.0 


0.383 ± 0.004 


13.8284 ± 0.0003 


-14.9005 


-0.689 ± 0.004 


0.69 



Table 1: Numerical results for the C F T R n F term of the R ratio. The three unphysical 
contributions, called 4-parton, 3-parton, and analytic, correspond to the three lines on the 
r.h.s. of eq. (3.30) respectively; their sums are reported under the column "total". The 
pulls are defined in eq. (4.65). 



increased, which implies that the integrand gives both positive and negative con- 
tributions. As is well known, when such contributions are close in absolute value, 
as when choosing y max = 0.6, the computation of the integral is affected by a rela- 
tively large error; in fact, the worst agreement with the exact result is obtained with 

2/max = 0.6. 

We finally computed a proper dijet total rate, by reconstructing the jets using 
the JADE algorithm with y cut = 0.1 (for a discussion on e + e~ jet algorithms, see 
e.g. ref. [40]). Using the same normalization as in ref. [19], we find for the NNLO 
contribution proportional to n F 

1.7998 ±0.0016 with y max = 0.6 , (4.67) 

which is in nice agreement with the result of ref. [19]. We point out that the small 
error in eq. (4.67), resulting from a run with lower statistics with respect to the 
results in table 1, is obtained with the value of y max that gives the worst convergence 
performance in the case of the total rate. We verified that, as for total rates, dijet 
cross sections are independent of y max within the error of the numerical computation 
(for example, with y max = 0.4 we obtain 1.7992 ± 0.0015). 

5. Comments 

We have shown in sect. 4 that our master subtraction formula, eq. (3.30), allows us 
to cancel analytically the soft and collinear divergences relevant to an NNLO com- 
putation. By using the ^-prescription, the subtraction counterterms rr_2, rf-i, and 
rv^ are constructed from the basic kernels that account for the matrix element sin- 
gularities. The successful implementation of the subtraction procedure also requires 
a sensible definition of the phase spaces used to integrate the counterterms. For con- 
venience, we collect here the phase space parametrizations we used, relabeling the 
partons where necessary, in such a way that those involved in the singular branch- 
ings have always the smallest labels. For the fr_ 2 counterterm to the double-real 
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contribution we used eq. (4.9) 



d*£ 2 = rf$ n (123) d$ 3 (l, 2, 3) , (5.1) 

relevant to topology 2 (topology 1(81 was not relevant to the case studied in this 
paper); we could also use eq. (4.12), which is fully equivalent and particularly suited 
to integrate kernels with a structure analogous to that of strongly-ordered limits. 
For the rr_i counterterm we used eq. (4.17) 

dKU = rf<Wl2) ^ (l " Jk) e *Ml,2) (5.2) 
= tten+idSj 1 . (5.3) 

Equations (5.2) and (5.3) give an explicit expression for d^&J 1 , which needed not be 
specified in eq. (3.28). Finally for the counterterm to the subtracted real- virtual 
contribution we used eq. (4.48) 

d$ n+1 = d$ n (12)— ^d$ 2 (l,2). (5.4) 

The factor that depends on s™ x in eq. (5.2) is crucial for the correct definition of the 
subtracted real-virtual contribution. The quantity s^ 3 * depends on the kinematics 
of the hard system whose phase space is <i$ n+ i(12). Such a dependence is in general 
non trivial: the difference in the results of eqs. (4.55) and (4.56) is due to it (notice 
that all of the 1/e poles are affected; thus, the failure to include the s^-dependent 
regulator in eq. (5.2) would prevent the cancellation of singularities). The dependence 
of s™ x on the kinematics of the hard system which factorizes needs not be known 
analytically for the computation of rr-idQ>~+ 2 an d for the definition of rv <s - ) . It must 
be known when defining fv^ s \ since this quantity has to be integrated analytically. 
However, in such a case the dependence is expected to be relatively simple (see 
eq. (4.49)). This is so because the subtraction term rv^ always corresponds to the 
kinematics of a strongly-ordered limit, in which the branching of partons 1 and 2 
will be followed by another two-parton branching (in the case studied in sect. 4, we 
had q gq 3 and g — > q^q'^ the case of topology 1<8>1 will clearly be even simpler). 

The phase spaces in eqs. (5.2) and (5.4) will serve to integrate the kernels. 
This implies that, for a given kernel, two different integrals will have to be computed, 
because of the presence of the s^^-dependent regulator in d$>~\. 2 - This can clearly 
be avoided by inserting the regulator also in eq. (5.4). It should however be stressed 
that, in doing so, we are forced to insert an analogous regulator (s^-dependent) 
in eq. (5.1). Thus, the advantage of having to compute half of the integrals rele- 
vant to the kernels may be lost because of the additional complications in the 
computations of the integrals of the kernels. 
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6. Conclusions 



In this paper, we have proposed a framework for the implementation of a subtraction 
formalism that allows us to cancel the soft and collinear singularities which arise in 
the intermediate steps of a perturbative computation at the next-to-next-to-leading 
order in QCD. The strategy is analogous to that adopted at the next-to-leading order, 
which is based on the definition of subtraction kernels whose form is both observable- 
and process-independent. We have introduced a systematic way (the ^-prescription) 
for defining the kernels, that results naturally in a two-step procedure. 

The first step is the definition of the subtraction counterterms for the double- 
real contribution. They are of two different types, denoted by fr_2 and rf_i; roughly 
speaking, the former collects all the pure NNLO-type singularities, whereas the lat- 
ter is singular in those soft, collinear, and soft-collinear configurations that are also 
relevant to NLO computations. The ^-prescription guarantees the absence of double 
counting. In the second step, the term rr_i is summed to the real-virtual contri- 
bution; this sum is free of explicit 1/e poles, but is still divergent, and a further 
subtraction needs be performed. This subtraction is completely analogous to that 
relevant to NLO computations. 

Our master formula, eq. (3.30), thus achieves the cancellation of the soft and 
collinear singularities essentially by two successive NLO-type subtractions; the first 
defines the subtracted real- virtual contribution, and the second removes its remaining 
singularities. The singularities that cannot be possibly obtained in this way are all 
contained in the subtraction counterterms fr_ 2 (the singularities of the double- virtual 
contribution are of no concern for the subtraction procedure - we assume the relevant 
matrix elements to be available). 

Equation (3.30) is, as it stands, not sufficient for an actual numerical compu- 
tation, since suitable choices of the phase spaces involved are necessary. We tested 
our subtraction formalism and the corresponding phase space choices in the context 
of a simple application, the calculation of the contribution proportional to C F T R 
of the dijet cross section in e + e~ collisions. Although this is clearly a simple ex- 
ample, it is, to the best of our knowledge, the first application in which NNLO 
process-independent subtraction counterterms have been constructed and integrated 
over the corresponding phase spaces to achieve an explicit cancellation of soft and 
collinear singularities. 

We point out that, although the phase-space parametrizations adopted in this 
paper are seen to induce reasonably simple analytic integrations, this is in part 
due to the relative simplicity of the kernels needed in the computation of the C F T R 
part of the e + e~ dijet cross section. More complicated kernels may require different 
parametrizations, which we did not consider in the present paper. In particular, the 
two NLO-type successive subtractions are liable to be simplified with respect to what 
done here. In any case, the general subtraction formula of eq. (3.30) will retain its 
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validity independently of the specific phase-space parametrizations adopted. 

There are a few difficulties that we did not address directly in this paper; they 
are not difficulties of principle, but may pose technical problems. One interesting 
feature will emerge when computing the C F T R part of the three-jet cross section in 
e + e~ collisions, namely the interplay between the 2 and 1<8>1 topologies which need 
be treated separately in the definition of the subtraction kernels. The case of colli- 
sions with one or two initial-state hadrons will also require a more involved notation 
(although the subtraction procedure will be basically unchanged). In general, it is 
clear that there is a fair amount of work to be done before the subtraction scheme 
proposed here is of any phenomenological use. This implies not only the definition, 
and integration over the phase spaces, of the universal subtraction kernels for all of 
the partonic branchings possibly occurring at NNLO, but also the computation of 
more two-loop amplitudes, which is a necessary condition for process-independent 
formalisms to be more convenient than observable-specific results. 
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A. The three-parton kernel K^h 

I1I2I3 

In this appendix we give more details on the construction of the three-parton kernel 
Kj\' q ■ According to eq. (4.5), the singularities due to the branching q — > qq'q' are 
described at NNLO by six terms. 

We start by showing that, of these six terms, only rr(CC) and r(C ®C) survive. 
Let us consider the limit in which the qq pair becomes soft. In this limit the matrix 
element squared behaves as [33] 

rp n 

rr(SS) = (47r«^) 2 — ^- ]T (M a3 T, • T, .H,,.,,J 

1^1 -K 2 ) ij=3 

In the limit in which q\ and q 2 become collinear the singular behaviour is given by 

rr(C) = ^a^)^P qq ,Mk T )T g Z... an 

= (-^ + Mi - *)^) r 9 z... an , (A.2) 

where T ga3 „ ,a n is the spin polarization tensor obtained by squaring the amplitude 
■Mga 3 ...a n and summing over colours and spins except the spin the gluon. 

Let us now take the double-soft limit of eq. (A.2). By using soft-gluon insertion 
rules we find 

rp n 

rr(CeSS) = -(^aX) 2 -^- £ <AW-d T * ' r ilM*-«J 

kf I , sk Ti jk TV \ kj 

x k t .(k 1+ k 2 ) \ 9fiu + 4z[i " z) ^r) k r (h+k 2 ) 

rp n 

= -(47^)2 £ (At a3 ... a J Ti • Tj\M a3 ...a n ) 

V 1 2/ i,j=3 

k^ k v 

" ( ^ r * 29 "" " 2KM M*7+*5) ■ (A ' 3) 

where we have used k\ = —2z(l — z)k±-k 2 . 

We now want to show that eq. (A. 3) is equivalent to eq. (A.l). We recall that 
the transverse vector k T is defined through the Sudakov parametrization 

fcf = Z pf> + k£-J?Z— n P, 
2zp-n 

k » =(1 - z) ^- K -- ( -^— n », (A.4) 
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where p is the collinear direction and n is an arbitrary additional vector used to 
define the collinear limit. Performing the replacement 

k T -> zki - (1 - z)k 2 (A. 5) 

we neglect terms proportional to the collinear direction, that vanish by gauge in- 
variance when contracted with the on shell matrix element. Gauge invariance, or 
equivalently, the fact that the eikonal current is conserved, can be further exploited 
to freely exchange k^ <-> —k 2fM in eq. (A. 5); we can thus perform the replacements 7 

kx/x fci/i k Tu > k 2v or 2k T ^k TU > k\^k 2v k\ v k 2ii (A. 6) 

thus showing the equivalence of eqs. (A.l) and (A. 3). We conclude that rr(C®SS) = 
rr(SS) and rr(C ®C®SS)= rr(CC © SS). 

We now turn to the computation of K~,%„ , which appears in eq. (4.22). As 
mentioned there, this quantity corresponds to the branchings q — > q%g, g q[q' 2 , 
which are described by the polarized Altarelli-Parisi kernels 



^(z,k T ) = ^zP%(±,k T ) , (A.7) 



By using the gluon polarization tensor 



-g^ + 4z(l - z) T T 



U2 

fVrp 



(A.8) 



d$(k) = -g» l ,+ k » n \ +Kn » (A.9) 



we can compute 



K S'l = — P£{*y + fc) dQ(h + k 2 ) p% ( k T ) d^(h + k 2 ) . (A.io) 



~ 1~ V-l -r -2, Ht, V m -r n, 2; , M - -t y u, v/3 

The vectors q T and fc T in eq. (A.IO) are the transverse vectors involved in the q q^g 
and g — > q[q' 2 collinear splitting respectively. By using the Sudakov parametrization 

tf = ^ + A£.-^l" < = 1,2,3, (A.ll) 

2p-n 

p being the collinear direction, q T and A; T can be identified with 

Qt — kx3 , kx = z 2 kxi — z\k T2 . (A. 12) 

We have 

k\ = -ziz 2 s 12 , q\ = -z 3 (zi + z 2 ) (si 3 + s 23 ) , (A.13) 
2g T -£; T = (^i + 22) (^2^13 - ^23) , (A.14) 



7 Note that the argument works in the same way starting from k T — > (1 — z)k\ — zk 2 - 
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where terms proportional to S12 in q\ and 2q T -k T have been neglected. Notice that 
d { ;J(h + h) = - 9lw + E^+P^t + (k Tt ) (A.15) 

and the 0(k Ti ) terms can be neglected in the evaluation of Kj^fa ■ Performing the 
algebra we obtain the result presented in eq. (4.25). 



B. Phase spaces 

This section collects the formulae relevant to the phase spaces used elsewhere in this 
paper. We refrain from giving the details of the derivations, which can all be worked 
out by using the techniques of ref. [35]. The phase space for the final-state particles 
of the process 

Q — > fci + ••• + *;„ (B.l) 

is denoted by 

d<S> n (Q; k u ...,k n ) = (2tt)^) (q - £ kij J] (27r)d -i2 fc o » ( R2 ) 

where d is the number of space-time dimensions; with abuse of notation, we may not 
list explicitly on the l.h.s. those momenta not involved in the necessary manipula- 
tions. The (d — l)-dimensional Euclidean measure is 

d d - l k t = \ki\ d - 2 d \ki\ dQ!f~ l) . (B.3) 

The angular measure is 

dftW = Q^- 1 ) ( S in ^_x) d " 2 dOa-! (B.4) 
= rt d ~V ( S i n Q^f- 2 dOj-! (sin e d - 2 ) d - 3 de d - 2 (B.5) 

= rt d -V [] (sin O^f- 1 - 1 d6 d ^ , (B.6) 



with 
and 



O<0i<2tt, 0<9i<n, % = 2, . . . , d - 1 , (B.7) 

s I *** = mry (a8) 

The recursive definitions of eqs. (B.4)-(B.6) can be continued analytically to complex 
values of d. Equation (B.l) can be rewritten as follows 



Q — ► q + k m +i H vk n 

> k\ + • • • + fc m , 
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where we have introduced the <i-vector q 



ki H h k n 



(B.10) 



Equation (B.9) suggests to decompose the n-body phase space in the following way 



d$ n (Q; k 1 ,...,k n ) = — d$ n ^ m+1 (Q; q, k m+1 , ...,k n ) d$ m (q; h, . . . , k m ) , (B.ll) 



which simply follows from inserting the identity 



(B.12) 



in eq. (B.2). In the context of higher-order QCD computations, the phase space 
decomposition of eq. (B.ll) matches the asymptotic behaviour of the matrix elements 
in the singular limits, which are expressed as a product of a kernel that collects all 
kinematic singularities, times a reduced non-divergent matrix element. In an NNLO 
computation, up to three parton momenta will enter the kernel. Thus, we shall use 
eq. (B.2) with m = 2,3. 

We start from the two-body phase space d& 2 (q; k 2 ). In order to give a 
frame-independent phase-space parametrization, we express the angular measure in 
eq. (B.3) in terms of invariants. In order to do this, two reference light-like four- 
vectors need be introduced, which define a system of axes; when azimuthal integra- 
tion is trivial, one of these vectors becomes irrelevant. After defining the reference 
vectors (which we denote by n and t, the latter relevant only to azimuthal integra- 
tion), the angles which enter eqs. (B.4)-(B.6) are written in terms of four-momenta 
through Gram determinants (see sect. 2 of ref. [35]): 



A„(pi,...,p n ) 



p\ 

Pi- Pi 



PrPi 
Pi 



PnPl PnPl 



Pl-Pn 
Pl-Pn 

Pi 



(B.13) 



We also define the variables 



Sll 



ki-n 
q-n ' 
2k r k 2 . 



1,2, 



(B.14) 
(B.15) 



After a somewhat lenghty algebra we arrive at 

& d -V A 2 (n,q)^A 4 (n,k u k 2 ,t)^ 



d<S> 2 (q;k 1: k 2 ) 



A 3 (n, h,k 2 ) 



d-4 
2 



8(27r) d ~ 2 (A 3 (n,A;i,A; 2 )A3(n,g,t))-*" V ^{n,q) J 
x 5 (1 — z\ — z 2 ) 5 (q 2 — S12 — k\ — k^j dz\dz 2 ds\ 2 ds\ t , (B.16) 
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where the second (third) fraction on the r.h.s. is due to the azimuthal (polar) part 
of the angular measure. By carrying out the azimuthal integration, we get 



d-A 
2 



4(2vr)^ 2 V A 2 (n,g) J 
x 5 (1 — z\ — z 2 ) 5 (q 2 — S12 — k\ — k^j dz\dz 2 ds\ 2 . (B.17) 



Using eqs. (B.14) and (B.15) we obtain 

A 3 (n, h, k 2 ) 



A 2 (n,q) 



= Z\Z 2 q — z\k 2 — z 2 k l . (B.18) 



We now proceed analogously for the three-body phase space d§ 3 (q; k±, k 2 , k 3 ). In this 
case, we shall not need to consider massive partons, and we thus set kf = from the 
beginning. Furthermore, one of the two azimuthal integrations will always be carried 
out trivially in our computations; thus, it is not necessary to introduce the auxiliary 
vector t, since its role can be played by one of the four-momenta ki (which was not 
possible in the two-body case). We get 



d$ 3 (q;ki,k 2 ,fa) - 



d-5 
2 



32(2tt) m - 3 V A 2 (n,q) J 
x 5 (1 - zi - z 2 - z 3 ) 5 (q 2 - S12 - si 3 - S23) dzxdz-zdzsdswdsisdsm , (B.19) 

where the variables Zi and have been defined similarly to eqs. (B.14) and (B.15). 
Upon explicit computation of the Gram determinants, we obtain 

A 4 (n, hi, k 2 , k 3 ) 1 / 

— r = -{2z 1 z 2 si 3 s 23 + 2ziz 3 s 12 s 23 + 2z 2 z 3 s u s 13 

A 2 (n,q) 4^ 

) . (B.20) 



_ 22_ 22_ 22 

— Z 1 S 23 ~ Z 2 S 13 — Z 3 S 12 



Equation (B.19) is fully symmetric in k i: and is therefore suited to treat the most 
general 1 — » 3 parton branching. There are cases in which such a symmetry is absent 
due to physical reasons (for example, a strongly-ordered double-collinear limit, i.e. 
C © C); in order to treat these, there are more convenient parametrizations than 
that of eq. (B.19). We use the factorization formulae eq. (B.ll) with n = 3, m = 2, 
and 

q — > p + k 3 

I (B.21) 

— >■ ki + k 2 , 

and we get 

d^ 3 (q;h,k 2 ,k 3 ) = ^-d$ 2 (q;k 3 ,p)d$ 2 (p;k 1 ,k 2 ), (B.22) 

Z7T 



38 



where for the phase spaces on the r.h.s. we can use eqs. (B.16) and (B.17). In 
particular, upon introducing the variables 



Zi = 



ki-n 




p-n 




212 




q-n ' 




q-n 


k\-n 


c 2 


k 2 -n 




p-n ' 


p-n 



S3(i2) = 2k 3 -p, (B.23) 



Ci = — , C 2 = — , (B.24) 
p-n p-n 

and using eq. (B.17) and eq. (B.16) we obtain 

fi(<*- 2 ) / A 3 (n,k 3 ,pY 



d-4 
2 



d$ 2 (q- k 3 ,p) 



4(27r) d - 2 V A 2 (n,q) J 
x 5 (1 - z 3 - z 12 ) 5 (q 2 - s 3 (i2) - p 2 ) dz 3 dz 12 ds 3{12) , (B.25) 



and 

Q(d-s) A 2 (n,p)^A 4 (n,A;i,A; 2 ,A;3)^ / A 3 (n,k u k 2 Y 



d-4 
2 



d$ 2 (p; k u k 2 ) 



^) d -\A 3 {nMM)^{nM + k 2) k 3 ))~ \ A2 ^) / 
x 8 (1 - Ci - C2) S (p 2 - sia) dddC 2 ds 12 ds 13 . (B.26) 



Upon explicit computation of the Gram determinants, we get 

= z 3 z 12 q 2 - z 3 p 2 , (B.27) 
C1C2P 2 , (B.28) 



\;{li.k:-,.j>) = . 2 : , 

A 2 (n,g) 
A 3 (n, fci,fc 2 ) 



A 2 (n,p) 

which are analogous to eq. (B.18). From eqs. (B.23) and (B.24) we also get 

z 12 = Zl + z 2 , Q = 7—^ — • (B.29) 

l — z 3 

We can introduce the azimuthal variable x through 

S13 = (S123 " si 2 )(l - C2) + j^C, (B.30) 

1-2:3 



S23 = (S123 - S12K2 - , Sl2 £ , (B.31) 

1-23 

where 8 

t = 2^3(2(1 - C2)(si 23 (l - z 3 ) - s 12 ) x + 2:3(2(2 - l)v^- (B.32) 

8 Note that the variable £ is directly related to the combination z 2 Si 3 — Z1S23 appearing in ti2,3- 
We have z 2 si 3 - z\s 2 3 = y/s^t,- 
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Using these definitions we rewrite the three-body phase space of eq. (B.22) as follows 



d-4 



x (C2 (1 — C2)) d;24 (1 - x 2 ) d 2 J ds l2 dz 3 d( 2 dx . 




) 



(B.33) 



We stress that eqs. (B.19) and (B.33) are identical, i.e. no approximation has been 
made in the latter. In fact, eqs. (B.33) could be obtained from eq. (B.19) by simply 
changing integration variables. 
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